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We investigate the QMCS in structure of the eigenfunctions, corresponding to mixed type 
classical dynamics in smooth potential of the surface quadrupole oscillations of a charged 
liquid drop. Regions of different regimes of classical motion are strictly separated in the 
configuration space, allowing direct observation of the correlations between the wave function 
structure and type of the classical motion by comparison of the parts of the eigenfunction, 
corresponding to different local minima. 



§1. Introduction 

The deformation potential 

U{ao,a2) = Y^ Cmn{al + 2a^)™a[f (6ai - a, 
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where ao and 02 are the internal coordinates of the drop surface 

R{e, ^) = i?o{l + ao>2,o(^, + a2[Y2,2{e, + Y2,-2{^, ^)]} 



(1-2) 




Fig. 1. The potential shape for W - 
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describes surface quadrupole oscilla- 
tions of a charged liquid drop of 
any nature, including atomic nucleP 
and metal clusters,^' containing specific 
character of the interaction only in the 
coefficients Cmn- Expanding to 
fourth order in deformation variables 
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V2C: 



„^ -ao,y= 02 (1-3) 

and assuming equality of masses for the 
two independent directions, we get the 
one-parametric Ca^-symmetric Hamil- 
tonian 



H 



2m 



+ Uix,y,W) (1-4) 



U{x,y,W) 



[x' + y^) + xy' 



-x' + {x' + y^f 



(1-5) 
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where W = 9CiJiCioC2o). 

The potential ()15|) is a generalization of the well-known Henon-Heiles potentiaP^ 
with one important difference: motion in (jl-Sf) is finite for all energies, assuring ex- 
istence of the stationary states in quantum case. For W > 16 the potential energy 
surface has seven critical points: four minima (one central and three peripheral) and 
three saddles. We consider in detail the case W = 18, when all four minima have 
the same depth E = and the saddle energies are Es = 1/20736. Critical energy 
of transition to chaos equals Es for the peripheral and roughly Es/2 for the central 
minimum, so we will be interested in the energy range Es/2 < E < Es, where the 
classical motion is chaotic in the central and purely regular in peripheral minima, 
resulting in what we shall call the mixed state.'^''^ 

We find numerically the spectrum En and the eigenfunctions ipn{x, y) for the Hamil- 
tonian Hl-4|) by the spectral method,"^*^ which implies the numerical solution of the 
time-dependent Schrodinger equation 



di^{x,y,t) 



^+[/(x,y,VF) 



^{x,y,t) (1-6) 



2m 

with the symmetrically split operator algorithiJ^ 

^(x, y, t^At) = e^^e-^^*^(^'^'^)e*^V'(a;, y, t) + 0{At^) (1-7) 

where exp {iAtV'^ /A'm)ip{x,y,t) is efficiently calculated using the fast Fourier trans- 
form. Initial wave function 

V'(x, y, t = 0) = ^ anipnix, y) (1-8) 

n 

is chosen to assure the convergence of 'ipi.x, y, t) in both the coordinate and reciprocal 
spaces and to avoid the degenerate states in the decomposition Hl-8|) . 
Having calculated ip{x, y, t) for < t < T, we obtain the spectrum En as the local 
maxima of 

T 

P{E) = i / dte''''Pit)wit) = Kl^^HE - En) (1-9) 



and the eigenfunctions ijjn{x,y) as 

T 

Mx, y) = J^J dt^x, y, t)w{t)e'^-' (MO) 



where 



P{t)= J dxdyr{x,y,t = Q)i:{x,y,t) (Ml) 

T 

5t{E -En) = ^j dtw{t)e''~^-^"^' (M2) 





and w{t) = 1 — cos(27rt/T) is the Hanning window function. 



QMCS in the Mixed State 
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§2. QMCS in the eigenfunction structure 

Quantum manifestations of classical stochasticity can be expected in the form of 
some peculiarities of concrete stationary stat^"^ or in the whole group of states 
close in energy)^~'2D Of course, it is not excepted that such alternative does not 
exist at all, i.e. the manifestations of the classical chaos can be observed both in the 
properties of separate states and in their sets. For example, comparing the eigen- 
functions, corresponding to energy levels below the critical energy of transition to 
chaos -Ec) with those corresponding to energy levels above Ec, drastic changes in the 
eigenfunction structure can be easily seen and analyzed. However, a question can 
arise if those changes are due to change in the type of classical motion or they are 
just a result of the quantum numbers change. Instead of comparing the eigenfunc- 
tions corresponding to different states, we propose to look for QMCS in the single 
quantum mechanical object — the eigenfunction of the mixed state. 
Cstj-symmetric Hamiltonian ()l-4p is invariant under rotation of 2ti/2> in the (x,y) 
plain and under reflection through x axis, so according to the transformation prop- 
erties 

'(/'(r, + 27r/3) = i?V(r, (2-1) 
'il){r,-ip) = a'4){r,ip) (2-2) 

the eigenfunctions can be divided into three types of different symmetry: Ai{R = 
1,0" = 1), A2{R = 1,0" = —1), and double degenerate E{R = exp (ib27r/3), o" = ±1). 
The initial wave function (jl-Bf) was chosen to excite only ^i-type states, and we 
computed the energy levels and the eigenfunctions for (|l-4() with W = IS and m = 
lO^'^, which corresponds to the main quantum numbers of order 10^ in the energy 
range Es/2 < E < Eg- Computations were made on a 1024 x 1024 grid of length 0.6 
(distance from the central minimum to saddles was 1/12 and 1/6 to the peripheral 
minima), number of time steps was 16384 with increment At = lO"^. 
Isolines of probability density and nodal curves for the eigenfunction corresponding 
to En = 0.444 X 10~^ = 0.92ii^5' are presented in Fig l2l31 As we did not compute the 
entire spectrum, we cannot point out the number n exactly, but we estimate it to 
be about 200. Correlations between the character of classical motion — developed 
chaos in the central minimum and pure regularity in the peripheral ones — can be 
easily seen comparing the parts of the eigenfunction corresponding to different local 
minima. 

§3. Conclusion 

We considered the quantum manifestations of classical stochasticity in the struc- 
ture of eigenfunction in the coordinate representation, corresponding to the mixed 
state in the potential of surface quadrupole oscillations of a charged liquid drop. 
Correlations between character of classical motion and structure of the eigenfunc- 
tion parts, corresponding to local minima with different type of classical motion, 
was observed in the shape of the probability density distribution and in the nodal 
curves structure, which gives direct and natural way to study the QMCS in smooth 
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Fig. 2. The eigenfunction shape Fig. 3. The nodal curves 



potentials. Numerical computations were made by the spectral method, which is a 
promising alternative to the matrix diagonalization method for the potentials with 
few local minima. 
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